Stickiness in a bouncer model: A slowing mechanism for Fermi acceleration 
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Some phase space transport properties for a conservative bouncer model are studied. The dy- 
namics of the model is described by using a two-dimensional measure preserving mapping for the 
variables velocity and time. The system is characterized by a control parameter e and experiences 
a transition from integrable (e = 0) to non integrable (e 0). For small values of e, the phase 
space shows a mixed structure where periodic islands, chaotic seas and invariant tori coexist. As 
the parameter e increases and reaches a critical value all invariant tori are destroyed and the 
chaotic sea spreads over the phase space leading the particle to diffuse in velocity and experience 
Fermi acceleration (unlimited energy growth). During the dynamics the particle can be temporarily 
trapped near periodic and stable regions. We use the finite time Lyapunov exponent to visualize 
this effect. The survival probability was used to obtain some of the transport properties in the phase 
space. For large e, the survival probability decays exponentially when it turns into a slower decay 
as the control parameter e is reduced. The slower decay is related to trapping dynamics, slowing 
the Fermi Acceleration, i.e., unbounded growth of the velocity 
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I. INTRODUCTION 

It is known that Hamiltonian systems are typical non- 
ergodic and non-integrable [ij . The phase space of such 
systems is divided into regions with regular and chaotic 
dynamics. These dynamical regions are connected by 
a layer, where regular or irregular motion, can or can- 
not mix, depending upon on the number of degrees of 
freedom of the system, as well properties of the limit- 
ing surface itself. Such a division leads to the sticki- 
ness phenomenon 0, which is manifested through the 
fact that a phase trajectory in a chaotic region passing 
near enough a Kolmogorov-Arnold-Moser (KAM) island, 
evolves there almost regularly during a time that may be 
very long. However, when an orbit resides in a chaotic re- 
gion far from the set of KAM regions, it moves chaotically 
in the sense that two nearby initial conditions apart from 
each other exponentially as the time evolves. Therefore 
the stickiness of phase trajectories has a crucial influence 
on the transport properties of Hamiltonian systems, and 
its relation to physical systems is one of the most im- 
portant open problems of nonlinear dynamics 0, [E] ■ Ap~ 
plications of stickiness can be found in astronomy [6j, 
fluid mechanics Levy flights [1| , also in biology Q , in 
plasma physics [lol . [T]| and many others. 

One of the main consequences of the influence of or- 
bits in sticky regime is observed in the transport prop- 
erties along the phase space. Therefore it may give rise 
to the following question: May sticky orbits influence 
the Fermi acceleration phenomenon? Fermi acceleration 
(FA) was introduced by the flrst time in 1949 by Enrico 



Fermi [12| as an attempt to explain the possible origin 
of the high energies of the cosmic rays. Fermi claimed 
that the charged cosmic particles could acquire energy 
from the moving magnetic fields present in the cosmos. 
His original idea generated a prototype model which ex- 
hibits unlimited energy growth and is called the bouncer 
model. The model consists of a free particle (making 
allusions to the cosmic particles) which is falling under 
influence of a constant gravitational field g (a mechanism 
to inject the particle back to the collision zone) and suf- 
fering collisions with a heavy and time-periodic moving 
wall (denoting the magnetic fields). The model is char- 
acterized by a control parameter e and has a transition 
from integrability e = to non integrability e ^ 0. A 
mixed structure of the phase space is observed for lower 
values of e and strong chaotic properties are present in 
the regime of large values of the parameter, say e > Cc 
where at Ec the system experiences a transition from lo- 
cal to globally chaotic regime (destruction of invariant 
spanning curves). 

In this paper we revisit the bouncer model seeking to 
understand and describe some transport properties along 
the phase space particularly focusing on the dynamics of 
sticky orbits. The model is described by a two dimen- 
sional, nonlinear and measure preserving mapping for the 
variables velocity of the particle and time at the collision 
with the moving wall. As the parameter e is increased, 
the number of islands in the phase space decreases. For 
the regime of high nonlinearity e ^ 1 , almost no islands 
are observed. The temporarily trapping dynamics due to 
the sticky regions are more often observed in the regime 
of small e where a mixed structure of the phase space 
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is present. We use the finite time Lyapunov exponent 
spectrum of the orbits and a statistical analysis of es- 
cape rates to investigate the influence of the stickiness 
in dynamics of an ensemble of non interacting particles. 
We therefore conclude that the stickiness present in the 
system acts as a slowing mechanism for FA. 

The paper is organized as follows: In Sec. |ll]the map- 
ping that describes the dynamics of the model is ob- 
tained. In Sec llin the numerical results are present which 
include the calculation of the finite time Lyapunov expo- 
nent and escape rates for the velocity as a function of e. 
Finally, the conclusions and final remarks are drawn in 
SecHvl 



THE MODEL, THE MAPPING AND 
CHAOTIC PROPERTIES 



II. 



moving wall. It leads to the following transcendental 
equation that must be solved numerically 

G(0c) = e cos((?!)„ -f- 0c) - e cos((?!)„) - Vn<t)c + ]^4>l ■ (2) 

If the particle leaves the collision zone case (ii) applies. 
The expressions are V* = —\JV^ + 2e(cos((/>„) — 1) and 
Ar„ = (t>u + 4'd+4'c with (j)u = Vn denoting the time spent 
by the particle in the upward direction up to reaching the 
null velocity, (p^ — \/Vn + 2e(cos(0„) — 1) corresponds 
to the time that the particle spends from the place where 
it had zero velocity up to the entrance of the collision zone 
at e. Finally the term (j)^ has to be obtained numerically 
from the equation F{4>c) = where 



. + 4'u + <t)d + 4'c)- 



-V* 
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We discuss in this section the procedures used to con- 
struct the mapping that describes the dynamics of the 
system. The model consists of a classical particle of mass 
m which is moving in the vertical direction under the in- 
fluence of a constant gravitational field g. It also suffers 
elastic collisions with a periodically moving wall whose 
position is given by y{t) — scos(wt), where w is the fre- 
quency and e is the amplitude of oscillation respectively. 

The dynamics of the system is made by the use of a 
two dimensional, nonlinear and measure preserving map- 
ping for the variables velocity of the particle v and time t 
immediately after a n*^^ collision of the particle with the 
moving wall. During the dynamics, two distinct kinds of 
collisions may be observed: (i) multiple collisions of the 
particle with the moving wall - those happening before 
the particle leaves the collision zone (the collision zone 
is defined as the region y G [—e,e]) - or; (ii) a single 
collision of the particle with the moving wall (causing 
the particle to leave the collision zone). Before writing 
the equations of the mapping, it is important to men- 
tion there are an excessive number of control parameters, 
3 in total, namely e, g and w. We may define dimen- 
sionless and more convenient variables as: Vn = Vnw/g, 
e = euP' I g and measure the time in terms of the number 
of oscillations of the moving wall 0„ = wtn ■ 

We assume that at the instant (j) € [0, 2tt] the position 
of the particle is yp{4>n) = ecos((/)„) with initial velocity 
> 0, which lead us to obtain the following expression 
for the mapping 
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K+1 = -V* +(f>c- 2esin((/)„+i) 
</)„+! = [(pn + ATn] mod(27r) 



(1) 



where the index c stands for the complete version of the 
model (the one which takes into account the movement 
of the moving wall) and the expressions for V* and Ar„ 
depend on what kind of collision happens. For case (i), 
i.e. the multiple collisions, the expressions are V* = Vn 
and AT„ = (j),. where (j)c is obtained from the condition 
that matches the same position for the particle and the 



The extended phase space for the whole version of the 
model considers four variables namely: {1) denoting 
the position of the moving wall; (2) Vp corresponding to 
the velocity of the particle; (3) Ep which is the energy of 
the particle and (4) the time t. The canonical pairs how- 
ever are: position and velocity (xw,Vp) and; energy and 
time {Ep,t). As the way the mapping was constructed, 
the variables used are not canonical ones therefore the 
determinant of the Jacobian matrix is 



Det J 



Vn 



Vn 
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which is clearly different from unity as it should be if 
the canonical pair was considered. However we may say 
that it preserves the following measure in the phase space 
d^i = [V + esiiY{(l))]dVd(l). 

A common version which is also present in the litera- 
ture is the so called simplified version. It was proposed 
many years ago [l^ as an attempt to keep the essence of 
the problem but at the same time allow numerical com- 
putations to be realized in a reasonable time when com- 
puters were far slower. Also it could reduce the complex- 
ity of the equations at a level that analytical calculations 
could be obtained. It assumes that the wall is fixed - 
so that the calculation of the time between collision does 
not evolve numerical solution of transcendental equations 
-, but at the instant of the collision, the particle suffers 
an exchange of energy and momentum as if the wall were 
moving. In this version, the extended phase does not 
consider more the position of the moving wall, because 
by definition it is fixed, causing the canonical pair to be 
the velocity and time. The mapping is then written as 



Vn+l = \Vn 
(pn+l = (pn - 



- 2fiSin{<j)n+l)\ 

2Vn mod(27r) 



(5) 



where the modulus function is introduced to avoid the 
particle to move beyond the wall. After a collision, if the 
particle has a negative velocity, we re-inject it back with 
the same velocity. For the simplified version and given 
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FIG. 1: Plot of the phase space for the bouncer model considering the control parameters: (a) and (e) e = 0.40; (b) and (f) 
e = 0.60; (c) and (g)e = 0.8; (d) and h) e = 1.20. 



the variables describing the dynamics are the canoni- 
cal pair, the determinant of the Jacobian is given by 
Det J = ±1. The simplified version of the model also al- 
low us to make a connection with the so called standard 
mapping. Defining J„ = 214,, K = 4e and 9n = (t>n+i + tt 
the simplified version is written as the standard mapping. 
The variation of the control parameter e leads the dy- 
namics to experience a transition from locally to globally 
chaotic dynamics as similarly observed in the standard 
mapping 14]. Indeed for e < Ec « 0.2429 the phase space 
has invariant spanning curves (also called invariant tori) 
and unlimited energy growth, which characterizes FA, is 
not observed. As the parameter e is increased, the fixed 
points become unstable and bifurcate for e > 1 {K > 4). 

The period- 1 fixed points are obtained solving the two 
equations simultaneously Vn+i = Vn = V* and (j)n+i = 
4>n — 4>* and are given by 

V* ^ nl , ; = 0,1,2,... (6) 

(j)* = arcsin (^^) , m = 0, 1, 2, . . . (7) 

Thus, there are windows of periodicity for the period one 



fixed points which depend on e. The linear stability for 
these fixed points are given by 

{2n f{p - 1)2 < 16e2 < {2n f{p - 1)^ + 4 , (8) 

where p = I — m. 

Figure [1] shows the structure of the phase space for the 
complete version of the bouncer as a function of the con- 
trol parameter e. The accuracy used to solve numerically 
both F and G was 10"^^ using the bisection method. As 
e is increased the stable regions (mainly marked by peri- 
odic fixed points) reduce leading the phase space to have 
large unstable regions. The regions of sticky are more of- 
ten observed for smaller values of e due to the existence 
of many islands in the phase space as compared to large 
values of e. Analyzing Fig. [T]we see that the phase space 
has a repeating structure in tt in the velocity axis. Thus, 
let us plot the phase space taking the morf(7r) for velocity. 
Such a plot is useful for observing the evolution of the 
fixed points and the possible trappings caused by sticky 
orbits. The control parameters used to construct Fig. [T] 
were: (a) and (e) e = 0.40; (b) and (f) e = 0.60; (c) and 
(g) e = 0.80; and (d) and (h) e = 1.20. For each figure 
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FIG. 2: (a) Plot of the FTLE for an initial condition chosen 
in the chaotic sea. (b) shows the evolution of the same initial 
condition of (a) for a plot of velocity against the number of 
collisions, (c) the zoom-in window of the previously selected 
area of (b) showing these trapping orbits in the phase space 
coordinates (V^, <^). 



a set of 100 different initial conditions were evolved in 
time until 10^ collisions with the moving wall. The ini- 
tial velocity was chosen such that its minimum value was 
higher than the stable region in y g [0, 2e]. 



III. NUMERICAL RESULTS 

This section is divided in two parts. In the first one we 
discuss the results for the Lyapunov exponent obtained 
at finite time while in the second we present our discus- 
sions and show results for orbits that survive longer the 
dynamics after being trapped by some sticky regions. 



A. Lyapunov exponents 

Let us start discussing our results for the positive Lya- 
punov exponent for chaotic components of the phase 
space. The Lyapunov exponent has been widely used 
to quantify the average expansion or contraction rate for 
a small volume of initial conditions. If the Lyapunov ex- 
ponent is positive, the orbit is said to be chaotic leading 
to an exponential separation of two nearby initial condi- 
tions. On the other hand, a non positive Lyapunov ex- 
ponent indicates regularity and the dynamics can be in 
principle periodic or quasi-periodic. The Lyapunov ex- 
ponents are defined as follows [11] (see for example [ll] 
for applications in higher dimensional systems): 





FIG. 3: (Color online). Plot of the FTLE distributions for 
several values of the parameter e. One sees two distinct peaks, 
a larger one representing the mean value of the Lyapunov 
exponent, and the secondary one, is due to orbits in stickiness 
regime. As e increases, the magnitude of the secondary peak 
decreases indicating that for higher values of e, less sticky 
orbits are observed. The control parameters used were: (a) 
e = 0.3; (b) e = 0.35; (c) e = 0.45; (d) e = 0.55; (e) e = 0.70; 
(f) £ = 0.80; (g) e = 1.0; (h) e = 1.2. 

where A", are the eigenvalues of the matrix M — 
Y\^^iJ{Vi,4n) and Ji is the Jacobian matrix evaluated 
over the orbit. 

In the dynamics of the bouncer model, chaotic and 
regular motion can coexist in the phase space, which in- 
troduces large variations and local instability along a ref- 
erence chaotic trajectory. Such variations, are related to 
alternations between different motions, in a qualitative 
way of saying, as well as chaotic and quasi-regular mo- 
tions. In order to characterize such peculiar variation 
dynamics, we used the Finite-Time Lyapunov Exponent 
(FTLE) [13. Once the trappings caused by orbits in 
stickiness regime happen just for a finite time, this tech- 
nique is useful to quantify the trapping effects. It was 
shown [17] that when the FTLE distributions present 
small values, it is related to existence of long-lived jets 
from a two-dimensional model for fluid mixing and trans- 
port. This can be understood, in a dynamic point of view 
as stickiness trajectories in the phase space. 

Figure shows the evolution of the FTLE, for an 
initial condition chosen in the chaotic sea, for e = 0.4. 
One sees a very irregular behavior along the time, al- 
ternating average contractions and repulsions, leading to 
and average value as A = 0.3078(1). In Fig. [Ifb) it is 
shown the evolution of the same initial condition of Fig. 
^Si) however plotted the velocity as a function of the 
number of collisions. It is clear in Fig. [Hb) the succes- 
sive trappings along the orbit, and how they "slow down" 
the energy growth, that characterizes the FA. Also, we set 
a zoom-in window in Fig. [^b) and plot the correspond- 
ing orbit in the phase space portrait {V,(j)), in order to 




FIG. 4: (Color online) Plot of the time evolution of initial conditions to reach a hole at Vhoic ~ 30. The control parameters 
used were: (a)e = 0.4, (b)e — 0.6, (c)e — 0.8 and (d) e = 1.2. Dark blue indicates long time evolution until reaching the hole 
while red indicates fast scape. White denotes the particle never escaping until 10^ collisions. 



identify some of these stickiness orbits in Fig. [H^c). 

To optimize the window of time to be used in the FTLE 
calculations, we have considered different lengths in sev- 
eral simulations. After some comparisons of the results 
we come up, based in fluctuations of the Lyapunov ex- 
ponents, to a finite time of 100 collisions that was then 
used to study the distribution of FTLE. It is known in 
the literature [l3|that the FTLE distribution has a Gaus- 
sian shape, where the large peak can be interpreted as 
the mean value of the Lyapunov exponent. If the system 
presents any periodic or quasi-periodic motion, besides 
chaos in its dynamics, the FTLE distribution can have 
a secondary peak in the region of very low value of the 
Lyapunov exponent. Such secondary peak is interpreted 
as sticky orbits along the dynamics evolution [l^, [l3] re- 
sponsible for trapping the dynamics. The distribution 
for several FTLE are shown in Fig. [3] for different con- 
trol parameter e as labeled in the figure. We can see 
from Fig. |3] that the secondary peak of the FTLE dis- 
tribution is more evident for small values of e. Just to 
have a glance of the influence of the second peak in the 



distribution represents up about 20% of the whole dis- 
tribution of Fig. Elb). The fraction of the distribution 
of the FTLE for the secondary peak decreases as e is 
increased. Such a result is expected because for higher 
values of e less islands in the phase space are observed as 
previously shown in Fig. [TJ 



B. Survival Probability and Escape Rates 

In this section we discuss results for orbits that survive 
until reaching a pre-defined velocity at which they are as- 
sumed to escape. To do that we consider the existence 
of a hole in the velocity coordinate of the phase space. If 
the particle reaches such a velocity or higher, its dynam- 
ics is stopped and a new initial condition is started. The 
introduction of the hole allow us to study transport prop- 
erties as well as characterize, through statistical analysis 
of survival probability and time-correlation decays, the 
influence of sticky orbits along the dynamics of the model 

MM- 
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FIG. 5: (Color online). Plot of the curves of Psurv for different control parameters. One sees the change of the behavior of 
Psurv as the parameter e is decreased. The control parameters used were: a)e = 1.4, b)e = 1.3, c)e = 1.0, d) e = 0.8, e)e = 0.6, 
f)e = 0.525, g)e = 0.475 and h)e = 0.425. 



To study the transport properties, we set a grid of ini- 
tial conditions equally distributed along the velocity and 
phase. Indeed a grid of 500 x 500 initial conditions with 
Vb € [e, 30] and 00 G [0, Stt] were considered. Then each 
initial condition was evolved in time up to the limit of 10^ 
collisions with the moving wall or until a hole placed in 
the velocity axis at T4ioic = 30 is reached. Figure |4] shows 
a plot of the initial conditions evolved until 10^ colli- 
sions with the moving wall or up to the particle reaching 
the hole. The color ranging from red (fast escape) to 
blue (long time dynamics) denotes the time (plotted in 
logarithmic scale) the particle spends until reaching the 
escape velocity. White regions denote that the particle 
never escaped. The control parameters used to construct 
the figures were: (a) e = 0.4; (b) e = 0.6; (c)e = 0.8 and; 
(d) e = 1.2. 

We see from Fig. SJa) where e = 0.40, that low ini- 
tial velocities spend large time accumulating energy un- 
til reach the hole at V — 30. Additionally one sees 
many stability islands where the orbits can get tempo- 
rally trapped and been released after a while. These tem- 



porally trappings are caused by sticky regions. Such dy- 
namical regimes can be visualized by the dark regions 
marked by blue color in Figs. |4l^b,c) whose control pa- 
rameters are respectively e — 0.6 and e — 0.8. When the 
control parameter e is raised, the particles reach the hole 
faster as we can see from Figs. |ll|b,c,d). In particular 
for Fig. \M,c) one sees that the first stability island dis- 
appeared. The stability regions are getting smaller and 
smaller as the control parameter e raises and from Fig. 
|4jd) they appear to be very small for e — 1.2. However 
even for a control parameter where the stability islands 
are small, we see that the sticky orbits are still present 
and indeed are marked by the dark blue color in the plot. 

The statistics of the cumulative recurrence time dis- 
tribution which is obtained from the integration of the 
frequency histogram distribution for the escape can also 
be obtained. To do that we consider now that the es- 
caping velocity is set as T4oio = 100 although any other 
velocity could be considered. Their cumulative recur- 
rence time distribution is also called survival probability 
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TABLE I: Exponents obtained from numerical fitting for the 
curves of Psurv for different values of e. 



and is obtained as 



1 ^ 



in), 



(10) 



where, the summation is taken along an ensemble of 
N = 10^ different initial conditions. The term Nrec{n) 
indicates the number of initial conditions that do not es- 
cape through the hole at T^ioie = 100 (i.e. recur), until 
a collision n. The ensemble of initial conditions was set 
for a constant velocity as Vb = 27r while 10® phase were 
distributed evenly in (/jq G [2.8,3.2]. 

It is known in the literature that if a system has fully 
chaotic behavior the curves of Psurv have an exponential 
decay [22i]. However, when a mixed dynamics is observed 
in the phase space, the curves of Psurv may present dif- 
ferent behaviors that may include: (i) a power law decay 
[2^ or; (ii) a stretched exponential decay [1^. For the 
bouncer model which has a mixed phase space the curves 
of Psurv may present either behaviors, depending on the 
parameter e and the set of initial conditions, as shown in 
Fig. [5l We see a transition in the behavior of the curves 
of Psurv as the parameter e is decreased. For large values 
of e as for example e = 1.4 and e — 1.3, the phase space 
has quite few islands and the chaotic sea is dominant over 
the dynamics. It is therefore expected an exponential de- 
cay in the curves of Psurv, as shown in Figs. [Ha,b). As 
the parameter e is getting smaller, more and more sta- 
bility islands appear in the phase space leading to the 
appearance of more and more sticky regions. With these 
stable regions around in the phase space, a change in the 
behavior of the curves of Psurv is expected. For values 
of e < 1 , we may observe a combination of decays in the 
curves of Psurv ■ Firstly the curves exhibit an exponential 
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FIG. 6: (Color online). Plot of — and —7 as a function of 



decay and suddenly they change to a slow decay that we 
observed to be described as a power law which marks the 
presence of orbits in stickiness regime 22 1 . 

Considering the curves of the survival probability 
shown in Fig. [SJ a numerical fitting can be made there- 
fore according to: (i) the exponential decay is given as 
Psurv{n) oc exp{nQ while; (ii) the power law decay is 
described by Psurv{n) oc where C and 7 are respec- 
tively the exponents for exponential and power law time 
decays. Table U shows the set of exponents for different 
values of the control parameter e. 

We see that as the parameter e decreases the exponen- 
tial decay of the curves of Psurv also suffer a change. The 
exponent C, decreases too as e decreases, a result which 
is quite expected given the periodic regions of the phase 
space are getting larger and larger. Figure [Sla) shows 
the behavior of the exponent ^ as a function of (e — ec). 
Looking at Fig. Hl^a) we see that the exponent C can be 
described by a power law of the type — C oc {e~ ecY and 
that the slope of the power law is given by z = 2.719(4). 

The exponent 7 however does not show the mathemat- 
ical beauty as observed for the exponent C,. The slower 
decay observed in the curves of the survival probability is 
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FIG. 7; (Color online). Plot of V as function of n for: (a) 
e = 0.8 and (b)e = 0.525. One can see two distinct growth 
exponents for Regular Fermi Acceleration and Sticky Fermi 
Acceleration. Such difference can be undesrtood as sticky 
orbits acting as a slowing mehanism for FA. 

indeed due to sticky regions present in the phase space. 
For our simulations, most of the slower decay was char- 
acterized as a power law. Indeed in the literature, it is 
known that the power law decay, for such cumulative re- 
currence time distribution for other dynamical systems 
[24 i25] which includes also billiards systems i26i - r2t| 
is set in a range of —7 G [1.5, 2.5] and that our results 
match this range. We stress however that the total un- 
derstanding and this behavior is still an open problem 
and extensive theoretical and numerical simulations, are 
required to describe its behavior properly. 

Let us now address specifically the assumption that 
stickness may affect the phenomenon of Fermi acceler- 
ation. Indeed the trapping dynamics of the particles 
around stable regions makes the unlimited energy growth 
slower than the usual. For a large set of initial conditions 
that lead the dynamics of the particle to present diffu- 
sion in the velocity, the average velocity V is described 
hy V (X \pn. However we expect the initial conditions 
that spend large time trapped in sticky regions lead the 
slope of growth to be smaller than 1/2. This is indeed 
true and figure [7] confirms this assumption. The curves 
shown in bullets in both Fig. [71[a,b) are named as Regular 
Fermi Acceleration (RFA) and were obtained for evolu- 
tion of the initial conditions which produce a fast decay 



in the survival probability (those along the exponential 
decay in Fig. [5]) and as expected, an exponent of = 0.5 
was obtained. On the other hand, the curves plotted as 
squares show the evolution of initial conditions chosen in 
the very final tail of the power law decay shown in Fig. [5] 
and are called as Sticky Fermi Acceleration (SFA). Power 
law fitting furnish slopes 0.398(7) for (a) and 0.400(1) for 
(b). These curves indeed give support for our claim that 
sticky regions slow down the Fermi acceleration. 
IV. FINAL REMARKS AND CONCLUSIONS 

The dynamics of the bouncer model was investigated 
by using a two dimensional measure preserving mapping 
controlled by a single control parameter e. For e = the 
system is integrable while it is non integrable for e 7^ 0. 
As soon as e increases, the periodic regions of the phase 
space reduce given rise to chaotic dynamics. Indeed for 
t> tc invariant tori are not observed in the phase space 
while periodic regions are still observed. The influence of 
sticky regions also reduces with the increase of e. Our nu- 
merical investigation of the FTLE spectrum distribution 
give support that trapping dynamics is often observed in 
the phase space and is confirmed by the secondary peaks 
of the FTLE distribution. The survival probability is 
characterized by two decaying regimes: (1) for strong 
chaotic dynamics, the decay is given by an exponential 
type while (2) it changes to a slower decay marked by a 
power law type when mixed dynamics is present in the 
phase space. Finally, according to the results shown in 
Fig 171 we see that when a strong regime of stickiness is 
present in the system, it acts as a slowing mechanism for 
FA. As with the survival probability, it would interest- 
ing to investigate whether the stickiness associated with 
mixed phase space in general models leads to a universal 
"slowing exponen" . 
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